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ABSTRACT.—The need to describe the relative intensity with which an animal uses different parts of its 
home range has been recognized for at least half a century. Such descriptions are particularly important 
for wide-ranging raptors with home ranges covering a variety of habitats. In studies of many taxa, the 
description of internal range structure is addressed by describing a core range of most intensive use. 
However, there is still no broadly accepted definition of a core or method of objectively estimating core 
ranges. Here, we propose that a core range can be usefully defined by the exclusion of excursive activity 
with the assumption that behavior differs between core and excursive activities. Two methods of ex¬ 
cluding excursive activity are presented for winter ranges of the Common Buzzard ( Buteo buteo) in 
lowland U.K. The first involves subjective exclusion of outlying locations, using the outermost discon¬ 
tinuity in the utilization distribution (UD). Incremental Cluster Polygons are used to produce the UD 
because this method provides the closest spatial relationship to the animal locations and the most clearly 
defined discontinuities. The potential for error or bias in this subjective method may often be unac¬ 
ceptable, particularly for home ranges which do not have well-defined core areas. The second method 
is a new application of incremental cluster analysis that objectively excludes excursive locations. The 
objective and subjective approaches are compared, and implications of core range definition in habitat 
and sociality analysis of raptors are explored in the context of published analyses on raptors and other 
taxa. 
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Estimation de rangos c\..n ale.: Una comparacion de tecnicas utilizando a Buteo buteo 

Resumen. —La necesidad de describir la intensidad relativa en que un animal utiliza las diferentes partes 
de su rango de hogar ha sido reconocida por lo menos desde hace medio siglo. Estas descripciones son 
particularmente importantes para las aves de presa que tiene rangos amplios con una gran variedad de 
habitats. En los estudios de muchos taxones, la descripcion de la estructura del rango interno es abor- 
dada mediante la descripcion del rango central como el mas utilizado. No obstante, no existe aun una 
definition aceptada del centro o de un metodo para estimar en forma objetiva los rangos centrales. 
Aqui proponemos que un rango central puede ser definido en forma util mediante la exclusion de las 
actividades de incursion. Dos metodos para excluir la actividad de incursion son presentados para los 
rangos de invierno de Buteo buteo en el Reino Unido. 

[Traduction de Cesar Marquez] 


Radiotelemetry has been used to study raptors 
for nearly three decades, providing data for many 
aspects of ecological research (Kenward 1985a). 
However, improvement in the collection and anal¬ 
ysis of data has been slower than the technical de¬ 
velopments in radio-tracking (Lance and Watson 
1980, Harris et al. 1990, Kenward 1992). The de¬ 


scription of the intensity with which animals use 
different parts of their home ranges presents a fun¬ 
damental problem (Hayne 1949). Animals gener¬ 
ally live in a spatially heterogeneous environment 
in terms of food availability, nest and roost sites, 
density of competitors, and other factors. There¬ 
fore, they tend to have one or more core areas of 
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intensive use in their home ranges (Kaufman 
1962) and it is likely that their behavior will differ 
in the core and outer areas of the home range. 
Many raptors are wide-ranging and their outer 
home-range boundaries may enclose habitat 
known to be avoided (e.g., Stahlecker and Smith 

1993) . Raptors may also make excursions from 
their usual range. Common buzzards (Buteo buteo), 
for instance, often make brief movements of up to 
20 km during their first year (Walls and Kenward 

1994) , before returning to ranges typically less 
than 1.1 km in diameter (Walls and Kenward 

1995) . Therefore, an analysis method for core 
ranges should exclude excursive locations (Burt 
1943) and areas within an outer home-range 
boundary that are avoided (White and Garrott 
1990). A number of methods have been proposed 
for estimating range cores, but none are widely ac¬ 
cepted. Therefore, comparisons between studies is 
generally not possible. 

The process of finding a core range can be split 
into three stages: (1) the description of the inter¬ 
nal range structure, giving nominal cores; (2) the 
choice of a percentage inclusion of radio locations 
that selects a biologically meaningful core for each 
individual; and (3) setting a standard core range 
(in terms of a standard percentage of locations) 
for a sample of animals to allow statistical compar¬ 
ison within the sample. Here, we present a com¬ 
parison of techniques used to describe the core 
ranges of raptors. 

Collection of Data 

Common Buzzards in southern Dorset, U.K. 
were instrumented with radio-tags at the nest just 
before fledging. Each radio-tag weighed 30 g and 
was mounted as a backpack with 6-mm-wide Teflon 
ribbon (Biotrack, 52 Furzebrook Rd., Wareham, 
U.K.). These tags had a life of up to 4 yr and a 
maximum range of 40 km from the ground and 
80 km from the air. To avoid disturbance of the 
study animals, buzzard locations were determined 
by triangulation from roadsides. Error associated 
with the locations was estimated at about 100 m. 
Standard 30 location home ranges (three locations 
per d for 10 d [Kenward 1987]) were recorded for 
122 buzzards from 1990-96. These included buz¬ 
zards aged between 1 and 4 yr. Data were collected 
in the nonbreeding season after the main dispersal 
period, when the buzzards had settled in relatively 
stable ranges. Data were analyzed using a modified 


version of RANGES V (Kenward and Hodder 

1996). 

Description of Range Structure 

Several different methods have been used for de¬ 
scribing internal range structure. The efficacy of 
these methods can be judged by their conforma¬ 
tion to the locations, including their ability to con¬ 
form to multinuclear cores. A further major re¬ 
quirement is efficiency. To give time for data col¬ 
lection on a sufficient number of animals, most 
projects require an analysis method that can esti¬ 
mate the range structure from a minimum number 
of locations per animal. The grid cell approach 
(Siniff and Tester 1965, Abies 1969, Voigt and Tin¬ 
line 1980, Samuel et al. 1985, Samuel and Green 

1988) conforms to location distribution. However, 
this method may require more than 150 locations 
to calculate a stable home range that does not in¬ 
crease in size as more locations are added (Don¬ 
caster and Macdonald 1991). The tessellation tech¬ 
nique proposed by Wray et al. (1992a) also re¬ 
quires a large number of locations. Contouring 
methods (Dixon and Chapman 1980, Worton 

1989) stabilize with fewer than 50 locations. How¬ 
ever, their accuracy of fit to the locations is influ¬ 
enced by dependency on an arbitrary grid and the 
use of parametric estimation functions (Spencer 
and Barrett 1984, Kenward 1992, Wray et al. 
1992b). Ellipses (Jennrich and Turner 1969) give 
stability with even fewer locations but conform 
poorly to the locations and can only provide one 
nucleus. Polygon-based techniques have stable out¬ 
er edges with 30 locations (Kenward 1982, Parish 
and Kruuk 1982, Kenward 1987); however, peeled 
convex polygons (Kenward 1985b, 1987) provide a 
poor fit to multinuclear or curved ranges (White 
and Garrott 1990). These problems may be avoid¬ 
ed with Incremental Cluster Polygons (ICP) (Ken¬ 
ward 1987). ICP analysis is based on forming 
groups of locations and separating outliers. Convex 
polygons drawn around the clusters provide a 
range outline that is not influenced by a grid or 
the position of outlying locations (Kenward 1987). 
Also, the outlines produced by elimination of out¬ 
lying locations stabilize at less than 50 locations 
(Kenward 1992). Therefore, we adopted ICP as the 
method for estimating internal range structure in 
this paper. 

Selection of a Biologically Meaningful Core Size 

The outline methods discussed can provide 
nominal core areas at any percentage inclusion of 
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locations. The second problem is to select a core 
that has biological significance. In the literature, 
the selection of core areas of most intense use has 
been largely subjective (Kenward 1985b, Harris et 
al. 1990, Wauters and Dhont 1992), or even arbi¬ 
trary (Mohr and Stumpf 1966, Anderson 1982, 
Wray et al. 1992b, Hohmann 1994). For instance, 
the core range has been defined as a 50% contour 
(Heikkila et al. 1996) or a 95% ICP core (Hulbert 
et al. 1996) but there is no biological basis for this. 
The subjective approach commonly uses the utili¬ 
zation distribution (UD) which is the polygon or 
contour area plotted against the percentage inclu¬ 
sion of locations (Van Winkle 1975). Identification 
of a discontinuity in this plot indicates the point 
where oudying fixes are excluded (Kenward 1985b, 
1987, Harris et al. 1990). ICP is particularly suitable 
for this method because it produces a stepped UD, 
unlike contour methods for which the plot tends 
to be smooth (Kenward 1987). 

An objective method for estimating a core range 
has been proposed by Samuel et al. (1985). How¬ 
ever, to achieve good conformity to the locations, 
their method depends on a large sample size of 
locations. In this paper, we analyze ICP ranges for 
buzzards to compare the subjective method (using 
the UD) with a new objective method for exclud¬ 
ing outlying (excursive) locations. 

Selection of Percentage Locations to be Included 
in Core Ranges 

Subjective Exclusion of Excursive Locations. 

When locations recorded for an animal include 
relatively long distance excursions (Fig. la), the 
outermost discontinuity on the slope of a utiliza¬ 
tion distribution can be very clearly defined (Fig. 
lb). In such cases, it is easy to visually select the 
point on the UD curve that indicates the size of 
the Excursion-Excluded Core (EEC). However, 
other individuals may have less well-defined cores 
(Fig. lc and Id). In this case, it is difficult to decide 
whether there is a core at 95% or 100% inclusion 
of locations (Fig. Id). EEC ranges were subjectively 
estimated from the UD curve for the ICP winter 
ranges of each of 122 buzzards we studied. 

Objective Exclusion of Excursive Locations. The 
separation of excursions from core ranges is based 
on the assumption that behavior differs between 
excursive and core activity. For example, buzzards 
and Northern Goshawks (Accipiter gentilis ) tend to 
fly for much of the time during excursions but 
make shorter, less frequent flights during more typ¬ 


ical foraging (Kenward 1977, Walls and Kenward 
1995). Therefore, if sufficient locations were re¬ 
corded, it is likely that the nearest neighbor (NN) 
distances between locations would form two fre¬ 
quency distributions representing core and excur¬ 
sive activity. In these hypothetical distributions, the 
mean distance between NN locations in the core 
would be expected to be smaller than the mean of 
NN distances where at least one of the locations is 
excursive. In most animal location data sets excur¬ 
sions are relatively rare. Therefore, the frequency 
distribution of NN distances would be expected to 
be negatively skewed with the positive tail indicat¬ 
ing excursive activity. 

The frequency distribution of NN distances was 
examined in a subsample of 10 buzzard ranges. For 
each of the ranges, the NN distance was calculated 
for all the locations (N = 30). To compensate for 
differences in the frequency distributions of NN 
distances between ranges, the distances were stan- 
d ardized by divi ding by h; for range; where h; = 
V (V^{a 2 xi + cr 2 yi }) and o- xi and cr yi are the standard 
deviations of the location coordinates for range; in 
the x and y directions (Worton 1989). As expected, 
the pooled frequency distribution of the standard¬ 
ized distances for all 10 ranges was highly negative¬ 
ly skewed, with many short distances and fewer lon¬ 
ger distances (Fig. 2). 

A transformation was sought that would high¬ 
light the excursive locations as outliers. There were 
only a small number of NN distances in each 
range; therefore, it was not possible to use tests of 
normality to seek an optimal transformation. In¬ 
stead, we sought the transformation that mini¬ 
mized the coefficient of variation (SD/ai) for the 
frequency distribution of NN distances. The effect 
of a number of logarithmic, reciprocal, and expo¬ 
nential transformations was examined on each of 
the 10 ranges. The transformation which most fre¬ 
quently minimized the coefficient of variation was 
generally the negative exponential transformation 
to q, where q = EXP (—0.5 (NN/h;) 2 ). This is the 
Gaussian kernel function of Worton (1989). This 
transformation was therefore applied to all 122 
buzzard ranges which were used to estimate cores 
with the subjective (UD) method. In order to sep¬ 
arate the excursive locations in each range, mini¬ 
mum excursion distances (MED) were estimated. 
These were the lower P = 0.05 (MED 0.05) and P 
= 0.01 (MED 0.01) percentiles of the normal dis¬ 
tribution fitted to the transformed NN distances. 
Clusters were then formed incrementally until the 



June 1998 


Core Range Estimates for Common Buzzards 


85 


a. c. 




95% 




Percent inclusion of fixes Percent inclusion of fixes 

Figure 1. (a) Range areas and fixes of a juvenile male Common Buzzard (Buteo buteo) (JM906). Hatched line indi¬ 

cates the 100% Minimum Convex Polygon, solid line shows the 90% Incremental Cluster Polygons, (b) Utilization 
distribution for incremental cluster analysis (ICP) of the winter range of male buzzard (JM906). There is a clearly 
defined discontinuity in the curve at 90% inclusion of fixes, (c) Range areas and fixes of a juvenile male Common 
Buzzard (JM939). Hatched line indicates the 100% Minimum Convex Polygon, solid line shows the 90% Incremental 
Cluster Polygons, (d) Utilization distribution for incremental cluster analysis (ICP) of the winter range of juvenile 
male buzzard (JM939). There is no clearly defined discontinuity in the curve. 


NN distance of the next fix to be added would have 
exceeded the MED. Any locations beyond this level 
were treated as excursive. The area of the convex 
polygon around the fixes in each cluster was 
summed to estimate the core range area. 

Core percentages and areas obtained with the 
two probability levels of the MED method were 
compared with those estimated by the subjective 
(UD) method. Results obtained by the UD and the 
MED methods were similar when MED = 0.05 and 
often larger when using MED = 0.01. For instance, 


buzzard JM939 had a 96% core at MED = 0.05, 
similar to the tentative core assigned with the UD 
method (Fig. Id) but had a 100% core at MED = 
0.01. Buzzard JM906 had a core at 90% inclusion 
of locations with the UD method (Fig. lb) and the 
MED method (0.05 and 0.01). Core areas of the 
122 buzzard ranges estimated by the UD method 
correlated well with core areas estimated using 
MED = 0.05 (r = 0.82, N = 122), and MED = 0.01 
(r = 0.89, N = 122). However, the distribution of 
the points around the 1:1 line showed that with 
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Distance of each fix from its nearest neighbor 
(NN) / hi (standardizing factor) 

Figure 2. Frequency distribution of nearest neighbor 
(NN) interfix distances pooled from locations (N = 296) 
from 10 winter ranges of Common Buzzards. NN dis¬ 
tances for each range are divided by a smoothing factor 
(h). Note that the y-scale is truncated. 


a. 



Core areas estimated using utilization 
distribution curves (ha) 


MED = 0.05 all except one range were close to the 
line, whereas with MED = 0.01 the majority of core 
areas were larger than those estimated subjectively 
(Fig. 3a and 3b). Therefore, MED = 0.05 was con¬ 
sidered appropriate for this sample of animals. 
This comparison showed that the MED method 
can give results which correlate well with the UD 
method, with the advantage that it is automated 
and does not involve subjective assessment for each 
animal. 

Setting a Standard Core Range Size for a Sample of 
Animals 

If a core is required with a standard percentage 
of locations for a sample of animals, it is advisable 
to select a percentage at which cores have been 
estimated for most of the ranges. We recommend 
setting a standard where 95% of the ranges have 
been cored, because this permits one range in 20 
to have many excursive locations without a dispro¬ 
portionate reduction in the sample core size. The 
resulting standard core percentage will only be 
larger than the core ranges of 5% of the sample. 
Therefore, few standard core ranges will include 
excursive locations. This is important because in¬ 
clusion of excursive locations greatly increases the 
core range area, and thus the variance of areas in 
the sample. A larger proportion of the cores may 
be larger than the standard percentage. However, 
this has a small effect on sample variance, because 


b. 



Core areas estimated using utilization 
distribution curves (ha) 


Figure 3. Comparison of core areas of 122 winter rang¬ 
es of Common Buzzards obtained using subjective (UD) 
and objective (MED) methods, (a) with MED P = 0.05 
and (b) with MED P = 0.01. The dotted line is the 1:1 
line. 


the removal of peripheral locations from these 
cores reduces the area much less than the removal 
of excursive locations. 

This process was applied to the sample of 122 
buzzard ranges. The number of ranges that had 
been cored increased as a function of the percent- 
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a. 



Percentage of fixes included in die 
Excursion Excluded Core (EEC) 


b. 



Percentage of fixes included in the 
Excursion Excluded Core (EEC) 

Figure 4. The core range size (percentage inclusion of 
locations) at which excursive activity is excluded by (a) 
subjective appraisal of the Incremental Cluster Polygon 
utilization distribution for buzzard home ranges (N = 
122) and (b) by the Minimum Excursion Distance (0.05) 
for buzzard home ranges (N = 122). 

age of core locations, for the objective (MED = 
0.05) and subjective (UD) coring methods (Fig. 4a 
and 4b). In each case, cores had been estimated 
for 95% of the buzzard ranges when 15% of the 
locations were removed (i.e., the excursion exclud¬ 
ed core contained 85% of the locations). With 
both methods, all the ranges reached cores that 
included at least 80% of the locations. 


Discussion 

In analyses of radio-tracking data, the choice of 
method for estimating home range size and struc¬ 
ture depends on the goal of the research (Voigt 
and Tinline 1980, Kenward and Walls 1994). Incre¬ 
mental Cluster Polygon analysis is particularly use¬ 
ful for identifying frequently used areas, as well as 
producing range structure statistics. Subjective 
(UD) and objective (MED) methods can be used 
to select the core ICP by excluding excursive lo¬ 
cations. However, the subjective choice of core 
from the UD may be difficult, particularly for rang¬ 
es that do not have a clearly defined core, and this 
has the potential to introduce error or even bias. 
For example, interpretation of the utilization dis¬ 
tribution could be influenced by prior knowledge 
of a typical percentage inclusion for cores in the 
sample of animals. In contrast the MED method 
provides a means to objectively plot a boundary 
that delimits the usual area of the study animal. 
The transformation applied to the frequency dis¬ 
tribution of NN distances and the choice of prob¬ 
ability level of the MED need to be tested with oth¬ 
er data sets. In the future, improvement in core 
range delineation might be achieved by plotting 
restricted edge polygons (Shekel 1954, Harvey and 
Barbour 1965, Voigt and Tinline 1980, Wolton 
1985) rather than convex polygons around the 
clusters. Further work is also desirable to test the 
efficacy of this approach for data other than the 
standard 30-location range, especially data includ¬ 
ing variable numbers of locations or collected dur¬ 
ing the breeding season. 

We suggest that excursion-excluded core ranges 
will provide important insights in behavioral ecol¬ 
ogy, especially in studies of sociality and habitat se¬ 
lection. In analyses of habitat use, a core range es¬ 
timator allows data to be viewed at three spatial 
levels: overall availability, familiar area, and the 
usual area sensu Burt (1943). The study area (usu¬ 
ally arbitrarily defined) can be used for overall 
availability. A range outline including all the loca¬ 
tions such as the Minimum Convex Polygon 
(MCP) (Mohr 1947) or a probabilistic contour 
(e.g., Worton 1989) can delineate the familiar area 
of the animal. Finally, a biologically meaningful 
core can reveal the usual area. The importance of 
including internal range structure in studies of 
habitat use has been demonstrated in analyses us¬ 
ing ICP cores. For example, in a study of Tawny 
Owls (Strix aluco ) in woodland patches, Redpath 
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(1995) showed that the owls had much larger MCP 
ranges where woodland was fragmented in com¬ 
parison to owls in continuous woodland. However, 
ICP cores (multinuclear polygons), were of a sim¬ 
ilar size in the different classes of woodland. 

The estimation of a core range may also be ex¬ 
tremely important in studies of behavioral inter¬ 
actions. Least overlap with conspecifics has been 
used to define core range boundaries (Auffenberg 
1978, Christian et al. 1986). However, a core range 
that can show if individuals regularly use the same 
areas, is more ecologically informative. For in¬ 
stance, when resources are very concentrated, in¬ 
dividuals may have overlapping home ranges, such 
as Northern Goshawks hunting near pheasant re¬ 
lease pens (Kenward and Walls 1994). Conspecifics 
that have overlapping outer ranges may show 
avoidance in their core (Kenward 1985b, Samuel 
et al. 1985, Harris et al. 1990). Intraspecific differ¬ 
ences in space use may also be masked by outer 
home range boundaries but revealed by the cores 
(Harris et al. 1990). 

ICP cores selected by the MED for the buzzard 
give a biologically useful estimate of the core 
range. We anticipate that the method presented 
here will prove useful for raptors and for other 
taxa. 
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